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Abstract 

The nucleation and evolution of bubbles are investigated in the model of an 0(3)-symmetric 
scalar field coupled to gravity in the high temperature limit. It is shown that, in addition to 
the well-known bubble of which the inside region is true vacuum, there exists another decay 
channel at high temperature which is described by a new solution such that a false vacuum 
region like a global monopole remains at the center of a bubble. The value of the Euclidean 
action of this bubble is higher than that of the ordinary bubble; however, the production rate 
of it can be considerable for a certain range of scalar potentials. 
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I. Introduction 



It is widely believed that the present universe was achieved through a number of phase 
transitions while the universe has expanded and cooled Q. A particularly interesting possi- 
bility is that of a first-order phase transition which is described by the formation and growth 
of bubbles . In the course of phase transitions in the early universe, two additional ingredi- 
ents should be considered: one is the gravitational effect [f| and the other is high temperature 
@, The inflationary scenario arose from the application of these ideas to cosmology is 
a synthetic attempt to answer the fundamental cosmological questions ||. These are the 
study of new phase bubbles in interiors of an old phase, but the opposite case, old phase 
remnants surrounded by the new phase, has also been an interesting subject 0. Recently 
inflation in the core of topological defects was analyzed in Ref. || ^ |10[ . Except for the 
last case, the order parameter of the model has been basically a real scalar field, thereby 
ignoring the difference between the continuous symmetry breaking and the discrete one at 
the time of bubble nucleation. Only for the evolution and collision of bubbles, this difference 
has been taken into proper account. However, the various species of solitonic defects, e.g., 1. 
topological or nontopological, 2. global or local, 3. cosmic strings, monopoles or textures, are 
determined by the continuous symmetry which the theory of interest holds [O]. Although 
soliton spectra involve complexities, the scenario of bubble nucleation is a simple one based 
on one real scalar order parameter, and continuous symmetry does not play a role for the 
formation of bubbles 0], whether it is due to quantum tunneling at zero temperature or 
due to thermal fluctuation at high temperature [|], || . Therefore it may be interesting to ask 
whether the structure of bubbles is affected by internal symmetry from the time of nucleation 
or not, and particularly how the internal symmetry coagulates a matter droplet inside the 
bubble. 

In this paper we consider a model of a scalar field with a global 0(3) symmetry and study 
the first-order phase transition at high temperature and in the presence of gravity. We first 
show that, in addition to the ordinary bubble of which the center is true vacuum, the model 
in a curved spacetime supports the 0(3) bubble that includes a matter droplet in its core as 



obtained in the same model in a flat spacetime 13]. The formation of the matter droplet at 



the core region of the new bubble is due to the winding between the internal space and that of 
spatial rotations, so it can be interpreted as a global monopole inside the bubble. Henceforth 
we will call it "monopole-bubble" in order to distinguish this bubble solution from the well- 
known ordinary bubble solution. This new bubble with a global monopole has two bubble 
walls: one is the bubble wall which distinguishes the false vacuum region outside the bubble 
from the true vacuum region inside it, and the other is the inner bubble wall which surrounds 
the false vacuum core of the global monopole. The long range tail of global monopole energy 
density is inversely proportional to the square of the circumference radius and affects the 
spacetime structure inside the bubble. Actually, it renders the spacetime between the inner 



and outer walls similar to the spacetime around the ordinary global monopole [|14]]: the 



spacetime has a repulsive nature, which indicates that the core does not collapse into a black 
hole even at the Planck scale, and the far region is a flat spacetime with a solid deficit angle. 
Although the total energy of the ordinary global monopole is infinite, in the present case it 
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is finite since the outer wall plays an important role of cutoff. 

The value of the Euclidean action needed to generate a bubble with a matter droplet is 
larger than that of an ordinary clean bubble; however, we observe that the decay rate for this 
monopole-bubble computed on the basis of the exponential formula can be quite considerable 
in comparison with that of the ordinary bubble for a certain range of the coefficients of the 
scalar potentials. Once bubbles are nucleated at high temperature in a curved spacetime, 
they start to expand because of some combination of processes when the environment keeps 
the temperature high [15] or the recovery of zero-temperature classical dynamics by the 
expansion of the background universe. We investigate the motion of the monopole-bubble 
by solving the coupled time-dependent field equations, neglecting the process of temperature 
changing. The outer bubble wall of the monopole-bubble grows similar to that of the ordinary 
bubble. The global monopole at the center of the bubble is stable when the phase transition 
scale is lower than the Planck scale. Moreover, this monopole-bubble also shares the same 
possibility of defect inflation at the Planck scale according to the arguments in Ref . 0, ^|, |1(| . 

The paper is organized as follows. In Sec. II we introduce the model and illustrate the 
nucleation of new 0(3) bubble solutions. In Sec. Ill we present the evolution of bubbles due 
to the gravitational force and briefly discuss the matter of defect inflation at the Planck scale 
at the core of the global monopole. Some concluding remarks are made in Sec. IV. In this 
paper we use the units c = h = k B = 1. 



II. Nucleation of Monopole-Bubbles 

Transition to the true vacuum state by quantum tunneling occurs through the nucleation 
of bubbles of the energetically favored phase. The nucleation rate per unit volume is T = 
Ae~ B , where B is evaluated by the Euclidean tunneling action and A is a prefactor which 
has units of energy to the fourth power. In this section we will consider a model with global 
0(3) symmetry and calculate the nucleation rate of bubbles at high temperature, specifically 
B of the model in a curved spacetime. 



A. Euclidean Solutions 



We begin with the action of an 0(3) symmetric scalar multiplet in the presence of Einstein 
gravity. Since the quantum statistics of bosons at finite temperature is formulated in the 
imaginary-time method which is described by the Euclidean theory with fields periodic in 
the Euclidean time with period (3 = 1/T, the Euclidean action at finite temperature is given 
by 

S = J Q P dt E J d^y/g^-j^R+^d^d^ + V^, (2.1) 



where a = a is an 0(3)-symmetric isovector (a = 1, 2, 3) with = ^fcjfcjf-. 
The equations of motion are 

1 d(-V) 
—d^gg^r) = — ^ (2.2) 
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= 87rGT M „ (2.3) 

where the energy-momentum tensor is 

= d,rd„r - {^-g pa d p ^d^ a + Vj . (2.4) 

Here we bring up the situation that the temperature is much larger than the inverse of the 
bubble radius, where the transition is described by time- independent fields. Both scales 
of symmetry breaking and temperature are usually assumed to be lower than the Planck 
scale Mpi = 1/yG; however, sometimes the supermassive scale when the transition scale is 
comparable with the Planck scale will also be considered when we examine the inner structure 
of bubbles in the very early universe. 

Although there is no rigorous proof in a curved spacetime when the 0(3) symmetric 
solutions with respect to spatial coordinates saturate the solutions with the lower Euclidean 
action, we choose to consider the bubbles which possess spherical symmetry. The form of the 
metric compatible with 0(3) symmetry can be written as 

ds 2 = (l - 2GM r {r) )e^dtl + (1 - 2GM r {r) ) "V + r 2 (d6 2 + sin 2 Ottf). (2.5) 
Using these coordinates the scalar field takes the following form, 

= (smn6 cosrrup, smnOsmnup, cos n8) <p(r), (2-6) 

where the spherical symmetry and the regularity at the origin allow only two cases: n — 0, 
and n — m — 1. 

Writing down the equations of motion by use of Eqs.( j2.5| ) and ( ^.6| ), we have three inde- 
pendent field equations: 

\ r ' dr 1 v r ' dr \ v r ' ) dr r l d<p 

ldS ^G(f), (2.8) 



r dr \dr 

1 dGM 



4nG 



r ' \dr J r 2. 



(2.9) 



r 2 dr 

where 5 n \ denotes the Kronecker delta. For the actual calculation, let us choose a sixth-order 
scalar potential such as 



A 



VU) = —U 2 + av 2 )((f) 2 -v 2 Y with < a < 1/2, (2.10) 

v 

which describes the potential in a broken phase. Although we choose a specific shape of the 
scalar potential, our argument in the following does not depend on the detailed form of the 
scalar potential and the existence of the new bubble solution with n — 1 is guaranteed under 



3 



any potential as long as it includes one false vacuum and one true vacuum. Here we only 
consider the transition from a symmetric vacuum to the broken vacuum, i.e. from de Sitter 
spacetime with the horizon H^ 1 = (8nGV(0)/3)~2 to Minkowski spacetime. 

The boundary conditions for nonsingular solutions of Eqs. (|2.7| ), ( |2.8|) and (|2.9| ) are 



— — (r = 0) = for n = or 0(r = 0) = for n = 1, 

dr 

4>(r -> H~ l ) = 0, and M(r = 0) = 0, 



(2.11) 



and we choose the normalization of t by setting 8{r — > H^ 1 ) = 0. The scalar field rapidly 
approaches the false vacuum for large r when the radius of the bubble is smaller than the de 
Sitter horizon H~ x . The above is our main interest; however, we will also comment on the 
cases in which the radius of the bubble is comparable to or larger than the de Sitter horizon. 
To find the behavior near the origin, we expand the variables by a power series, finding 



5{r) 
GM(r) 



ldV_ 
6 ~d~4> 



r 2 , (fiesc = 0(0) 



5 esc + t;G( 



9 \d(f) 
r ^n>.w)r :1 



) r 4 , 5 CSC = 5(0) 



(2.12) 



for the n = configuration, and 

<j)(r) 

5{r) 
GM(r) 



1 



iA(2 

2 2 
r 



a: 



5 + 2nG4> z r 
2nG(j)lr 3 



= 5(0) 



^(O) 
dr 



(2.13) 



for the n = 1 configuration. The constants (0 eS c, <5 eS c) and ((/>o, 5q) are determined by the 
proper behavior of the fields for large r. If the transition scale ran = ^4A(3 + 2a)v is 
smaller than the Planck scale and the size of the bubble is smaller than that of the de Sitter 
horizon H , we may meet no coordinate singularity inside the bubble. Therefore, the scalar 
amplitude <p(r) and a metric function 5(r) approach their boundary values exponentially in 
their asymptotic region, and Eq.(|2.9|) says 



GM(r) 



-nGXav r 



4^3 



(2.14) 



since the first two terms in the right-hand side of Eq. ( |2.9| ) are negligible and the last cosmolog- 
ical constant term is dominant at the false vacuum region outside the bubble. Figure 1 shows 
that 5(r) and GM(r) are monotonically increasing functions of r; however, the behavior of 
the scalar field <f> near the origin depends on whether n — or n — 1. 
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Figure 2: T\ profiles for fixed A = 1 and v/Mpi = 0.1. (a) thick-wall bubble of a = 0.3, (b) thin-wall 
bubble of a = 0.01. The dotted and solid lines correspond to an n = bubble and an n = 1 bubble, 
respectively. 
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As shown in Fig. 2, the n = solution is the well-known vacuum bubble solution where 
the energy is accumulated at the bubble wall. Although there remains some matter inside 
the thick- wall bubble (see the dotted line in Fig. 2-(b)), contrary to the case of the thin- wall 
bubble (see the dotted line in Fig. 2-(a)), such matter does not form an aggregate inside the 
bubble. The energy density of n — 1 bubbles, 



(2.15) 



is expressed by solid lines in Fig. 2. We easily read that a matter droplet is formed at 
the center of the n = 1 bubble due to the nontrivial local mapping between internal 0(3) 
symmetry and spatial 0(3) symmetry. From the boundary conditions of the scalar field in 

we can interpret this matter droplet global 



Eq.( |2.11|) together with the ansatz ( ^6| 
monopole. In a flat spacetime of dimension more than two, there is a no-go theorem that 
says the scalar fields described by the standard relativistic form of the Lagrangian do not 
support non-trivial static soliton solutions of finite energy [16]. Therefore, when we consider 
global vortices or global monopoles in the presence of gravity [JD], [TJ]], the introduction of 
a cutoff scale, for example, the horizon length, provides us a way to control the divergent 
quantities fl8 |. However, the global monopole created inside the n = 1 bubble is a finite 
energy static configuration since the long-range tail of the global monopole is tamed by the 
outer bubble wall. 
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The features described above do not depend on whether the bubbles lie in a flat spacetime 

or a curved spacetime; however, other characteristics due to gravity are also worth noting. 

We can read the gravitational effects on the shapes of bubbles in the weak gravity limit 

as follows. If we think of r time and 6 position of a hypothetical particle in a 

one-dimensional motion, the scalar field equation ( |2.7| ) can be interpreted as the Newton 

/ 2GM\ 

equation for a particle of variable mass (1 J subject to several forces: the first one is 

the conservative force from — V(<p), the second one a friction with time-dependent coefficient 

/ 2GMn d , / ox,, 2GM ,\ , , , , P , . 

(1 — ln(r e (1 ) , and the last one only tor n — 1 is a time- dependent 

v r s dr ^ r ' 

repulsion whose strength is extremely large at r = 0. Although it is difficult to prove 

analytically the existence of solutions in the presence of gravity, we can understand how 

and what kind of bubble solutions are supported by use of the terminology of Newtonian 

mechanics as was done in the flat spacetime case ||19|| . Let us consider the returning motion 

of a hypothetical particle both for n = and n = 1 solutions, i.e., from 0(r = 0) = e sc for 

the n = solution (or 0(r = r turn ) = tU rn for the n = 1 solution) to <p{r = oo) = 0. Since 

2GM 

GM(r) is an increasing function of r, the variable mass of a hypothetical particle (l j 

decreases as time r elapses. Furthermore, the gravitational effect due to the sum of S(r) and 
GM{r) decreases the time-dependent coefficient at the starting point for the n = solution 
and at the turning point for the n = 1 solution. Both effects involve the decrement of energy 
gained during the rolling of a hypothetical particle down to the minimum point of the effective 
potential — V(0), hence CSC (or (fitum) should become small after the inclusion of gravity (see 
Fig. 3). For an n = bubble, it makes the radius of the bubble smaller. For the behavior 
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of the n — 1 bubble, the hypothetical particle reaches the returning point, 0tum, at a time 
earlier than in the flat case. In order to reach the smaller turning point 0tum i n spite of the 
larger acceleration for a given 0o and for small r as shown in Eq. (|2.13| ), the initial velocity 
0o of the particle for n = 1 bubbles should be smaller than the value in a flat spacetime. 
Hence, the radius of the bubble becomes small, although the size of the matter droplet is 
increased by the gravity. The contraction of the outer walls of both n = 0, 1 bubbles can 
be understood through the attractive nature of gravity; however, the explanation for the 
larger matter droplet will be given later by mentioning the repulsive nature of gravity on the 
matter core. Looking back at Fig. 2, we can compare the radius of an n = bubble and of 
an n = 1 bubble. In a flat spacetime the radius of an n = bubble is always larger than 
that of an n = 1 bubble [|13j , and it is also true for a thick- wall bubble in a curved spacetime 



(see Fig. 2-(b)). For a thin-wall bubble in a curved spacetime where the global monopole 
stretches a long-range tail; however, the radius of an n = 1 bubble becomes smaller than that 
of an n = bubble. This fact will be justified by our analytic discussion under the thin-wall 
approximation in Sec. III. 

For an n = 1 bubble, let us look at the spacetime structure in the neighborhood of the 



point r = r turn where takes the maximum value 0t U rm i-e., — 

dr 



0. Equations (|2.8|) 



-^tur: 



and ( |2.9| ) around r = r turn = have approximate solutions such as 

5(r) « (Wrn (2.16) 
GM(r) » GM tum + 4vrG0 t 2 um ^ + ^GV r (0 turn )r 3 . (2.17) 
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The constants 5 tU m and M turn are fixed by 



-4nG I" '*-r(^) 2 (2-18) 



GM tura „ 4^fV^{I(l-^)(^)%fl^ + ,K-K(^,)}.,2.19) 

The above integrals are estimated as <5 turn ~ 47rG0^ urn and |GM turn | ~ |Gm#|. 

If we consider an n = 1 bubble with a thin wall when the phase transition scale m# = 



4A(3 + 2a)v is much smaller than the Planck scale Mpi, the first and the third terms in 
Eq. (|2.17|) can be neglected since V(4> ivaiL ) ~ and ur turn 3> 1- Substituting the above results 
into Eq. (|2.5D and rescaling the variables t and r as 

t _> (l-87rG0 t 2 ur n)^e- 5 -' n t (2.20) 

r _» (l-87rG0 t 2 ur „)^, (2-21) 

we obtain a metric after a Wick rotation, which describes the region inside the outer wall 
but outside the global monopole: 

ds 2 = -dt 2 + dr 2 + r 2 (l - 87rG0 2 urn )(d# 2 + sin 2 9 dip 2 ). (2.22) 

Although the actual metric is not completely flat due to the additional small terms in 
Eq. (|2.17|) , the observer inside the n = 1 bubble feels no gravitational force exerted by the 
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global monopole apart from the tiny effects from the monopole core and the energy dif- 
ference between the true vacuum v and the maximum value of the n = 1 bubble 4>tum, 
^(0tum) — V(v). This phenomenon can be explained by a Newtonian gravitational potential. 
The radial component of tension — T r also has a long range term such as 



r 



which cancels the energy density T\ in fl2.15| ) in the Newtonian limit of the Einstein equations: 
V 2 $ = 8ttG{T\ — T r r )^0 at r « r turn . However, since the metric (|2.22| ) describes a space with 



a deficit solid angle, if we consider a light signal propagating from a source to an observer, 
the observer inside the n = 1 bubble must notice the light bending due to the deficit solid 
angle A = 8irGv 2 , as is the case of a straight cosmic string. Thus a rough evaluation 
gives the angular separation 5(f ~ 8tcGv 2 ~few arcsec at a typical grand unification scale 
v ~ 10 w GeV, which can be observable. 

For the phase transition in the supermassive scale, the absolute value of M turn in Eq.( |2.19| ) 
is of the order of the Planck scale, |M turn | ~ Mpi, and then the first term in Eq. fl2.17| ) becomes 
considerably large. Let us discuss the structure of the spacetime in this case. Assuming the 
outer wall is extremely thin, i. e., the third term in the right hand side of Eq. ( |2.17| ) is negligible 
inside the bubble, we obtain a metric around <f)(r) ~ 0tum as 

rfs 2 = _(i_ S7rGv 2_ 2GM ^ yS^ dt 2 + (l_ g7rGv 2_ 2G?Mturll )- 1 t j r 2 + r .2 (d g2 +sill 2 ^2) 

(2.24) 

Here we make a crude assumption that the region between the inner wall and the outer wall 



(R m < t < Rn=x, see Fig. 2) is described by the above metric in Eq. ( |2.24| ), and the region 



inside the inner wall (r < R m ) is approximated to be de Sitter spacetime: 

ds 2 = -(l - ^7rGV(0)r 2 )e 25o dt 2 + (l - jjvrGV (0)^)^ dr 2 + r 2 (d6 2 + sin 2 ^ 2 ). (2.25) 

Equation ( p.8| ) says that the change of 5(r) can be neglected inside the monopole, which is 
supported by Fig. 1(b), and thus we take 5 ~ 5 turn under our assumption. From the conti- 
nuity of the metric and its first derivative with respect to r, we estimate the size of the global 



monopole (or equivalently the position of the inside bubble wall) as R m = l/w4A(3 + 2a)v 



and M turn = — y4A(3 + 2a)v < 0. It is confirmed by envisaging the integral in (|2.19| ): if we 
divide the integration domain into (0,R m ) and (R m , r turn ), and then substitute the values of 
the scalar amplitude, i.e., = in (0,R m ) and = tU rn in (R m , r tum), the integral of the 
core region has only the negative contribution to GM tmn . This result is consistent with the 
known result that the global monopole inside the n = 1 bubble does not form a black hole 
even at the Planck scale [Q. It can also be checked by the radial motion of a test particle 
governed by the geodesic equation 

d 2 r _ d /GM^ ^ 2GM^d5 
dr 2 dr \ 
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where r is the proper time of a test particle. Equation (|2.26|) , together with fl2.13|) , tells us 
that the acceleration increases in proportion to the radius around the center of the global 
monopole. Although it is a weak gravity case (v/Mpi = 0.1), Fig. 4 shows a typical example 
consistent with the above argument. Here a question arises: what is the structure of a 
spacetime manifold which is formed when the deficit solid angle is equal to or greater than 
47r? In the case of local cosmic strings, when the deficit angle is equal to or greater than 



2tt, a possible two-dimensional spatial manifold is described by cylinder or two sphere [pO 
However, it is an open question for the global monopole. 




Figure 4: The acceleration of a test particle d 2 r/d,T 2 versus the radius r for A = 1 and a = 0.1. 
The acceleration denoted by the solid line is always positive and increases at the core of the global 
monopole for a weak gravity case (v/Mpi = 0.1). 



B. Nucleation Rate 

We now turn to the evaluation of the nucleation rate of bubbles; specifically the values 
of the Euclidean action at two stationary points of n = (Bo) and of n = 1 (B\) bubbles 
are of our interest. Suppose that there exists a barrier between a local minimum of B Q and 
Bi, and the maximum value of B between B and B\ can be set to be -Btop- When each 
valley is so deep that the height of the barrier at the hilltop B top evaluated from the ordinary 
bubble point is larger than the difference of two local minima at the bottoms of valleys, i.e., 
Btop — B\ ^> B\ — Bq, then the total decay rate from a metastable phase into a stable phase 
per unit volume can be estimated to be the sum of the nucleation rate of each bubble 

r = r (o) +r (i). (2.27) 

It means that both n = and n = 1 solutions give distinct decay channels where each 
solution describes the nucleation of bubbles with the critical size. Moreover, if the tunneling 
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action for each bubble is larger than unity, the nucleation rate for the ra-th bubble takes the 
exponential form 



(2.28) 



On the other hand, if the tunneling action is of the order of unity or smaller, which corre- 
sponds to thick- wall bubbles in the high temperature limit, the exponential formula ( 2.28| ) 
is no longer valid. Here we take a heuristic viewpoint and keep our analysis on the basis of 
the above formula in Eq.( |2.28| ), notwithstanding the above possibility. When the first-order 
phase transition is considered in a curved spacetime, the background spacetime itself expe- 
riences time evolution, e.g., inflation in the metastable vacuum region and this expansion 
of spatial volume induces a sudden drop in temperature. Therefore, the shapes of bubbles 
or equivalently the configurations of the scalar field change, and hence the formula (|2.28| ) 
should not be applied. Here we discuss the probability to nucleate bubbles when the system 



is initially in the metastable phase and the temperature change can be neglected pl| . Under 
this restriction B n can be approximated to be the values of action fl2.1|) for a given n = or 
n = 1 bubble since we consider only the decay from a spacetime with a positive cosmological 
constant, V(0) > 0, into a spacetime with zero cosmological constant, V(v) = 0. 

We plot a dimensionless value B' n = (T/v)B n for various v/Mpi and a in Figs. 5 and 
6, respectively. Figure 5 shows that B n becomes small for large v/Mpi, which implies that 
the materialization of both n = and n = 1 bubbles is more likely at higher energy scales. 
This result is consistent with the fact that the gravitational contribution to the total energy 
is always negative, as we have already seen in Fig. 3. 
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Figure 5: The values of (T/v)B n , shown versus the vacuum expectation value of the scalar field in 
unit of the Planck scale, v/Mpi. The solid and dotted lines correspond to an n = bubble and an 
n = 1 bubble, respectively. 
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Figure 6: The values of (T/v)B n , shown versus the parameter a which governs change of the shape 
of bubbles. The solid and dotted lines correspond to an n = bubble and an n = 1 bubble, 
respectively. 

When we have a model of a first-order phase transition with several decay channels, an 
interesting question may be which decay channel is dominant, which is illustrated by the 
ratio of the two decay rates, r^/r ■ The prefactors A n usually are assumed to be of the 
order of m 3 H T at high temperature. If we neglect the difference between A and A\, then the 
relative decay rate is determined by the exponential factor, 



r(o) 



cxp 



—A 



(2.29) 



As expected, the n = 1 bubble solutions take higher values of action, i.e., B\ > B , regardless 
of the shape of the scalar potential (see Fig. 6) and of the strength of gravitation (see Fig. 
5). From Fig. 6, we read that B x — B' is large in the thin- wall limit (small a), because the 
leading contribution of action difference (v/T)(B[ — B' ) can be understood as the energy to 
support the global monopole in an n = 1 bubble, of which the long-range tail of the global 
monopole, T\ ~ v 2 /r 2 , consumes the energy proportional to the radius of the n — 1 bubble. 
On the other hand, if the bubble wall is relatively thick, B[ — B' becomes small, and then 
the nucleation of n = 1 bubbles is not negligible. When the relative ratio of the exponentials 
( |2.29[) is not far from the order of unity, we should take into account the ratio of prefactors 
Ai/Aq in order to determine the dominant decay channel. It is difficult to compute A n even 
for an n = bubble in a flat spacetime. In a curved spacetime it has units of energy to the 
fourth power and is expected to be of the order of vn? H T . To get some information on the 
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ratio of the prefactors, let us consider bubbles in a flat spacetime ||, [T3| 



(2.30) 

where the system for the fluctuations around each classical solution includes six zero modes 
(three from spatial translations and another three from spatial rotations). As explained 
before, Bi/B — 1 represents the ratio of energy to make a global monopole to that to 
generate a bubble, thus Bi/Bq tends to one in the thin- wall limit and a few in the thick- wall 
limit. If we just replace the values of the action in a flat spacetime to those in a curved 
spacetime, we obtain several values of /T^°\ as displayed in Table 1. 
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5.23xl0- 13 


1.47xl0~ 25 



Table 1. Values of rW/T (0) for A = 1 and vjM Pl = 0.1. 



We find that monopole-bubbles become more likely to be nucleated at high temperature and 
in the relatively thick-wall case. (Remember that larger a corresponds to a smaller potential 
barrier, which creates a bubble with a thicker wall.) Although B\ is always larger than B , 
there may exist some parameter region of the scalar potential where n = 1 bubbles cannot 
be neglected. 



III. Evolution of Monopole-Bubbles 

When we consider a first-order phase transition in the framework of a finite-temperature 
field theory with imaginary time, high-temperature bubbles are given by static solutions of 
the Euclidean equations, and they are also static solutions of the Lorentzian equations. It 
is obvious that the bubbles start to evolve immediately after their nucleation, so a way of 
description is to borrow the physics of combustion processes when the environment keeps 
the temperature high enough ||15|| , which is indeed the case in a flat spacetime. Once the 



gravity is taken into account in the early universe, the background universe is expanding and 
then it is rapidly cooled down to zero temperature. The motion of bubbles eventually follows 
the zero-temperature classical dynamics. An accurate bubble dynamics may be as follows. 
When the bubble larger than the critical size is nucleated, the detonation process induces 
the growth of it and simultaneously the region outside the bubble begins to expand due to 
the gravitational effect. In the next step the variation of temperature requires the inclusion 
of whole complex ingredients into the evolution procedure of bubbles; for instance, the tem- 
perature dependence of the the combustion processes, the effect of gravity on the classical 
evolution of bubbles, the change of the effective potential due to temperature and quantum 
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corrections, possible reheating, the generation of Goldstone bosons and so on. However, we 
already know what actually happens to n = bubbles: if the expansion rate of the universe 
is large enough, the motion of an n = bubble turns out to be that of a bubble governed by 



the classical equations of motion at zero temperature [|2|, |22|, Here we suppose that the 
above simplification for n = bubbles is applied to the case of n — 1 bubbles in a similar 
manner and concentrate our interest only on the evolution of n = 1 bubbles due to classical 
effects. 

Our task is now reduced to solving the Einstein equations and the scalar field equation 
by use of numerical analysis. In order to examine the evolution of n — 1 bubbles after the 
temperature decreases to zero, we should prepare the initial configurations. However, they 
should be different from our static n = 1 bubble solutions of the Euclidean equations because 
they are also static solutions of equations of the motion after a Wick rotation. The effect 
due to the change of temperature must be reflected when we prepare the initial conditions. 
Suppose that the various effects mentioned above give rise to the evolution of bubbles and then 
the structure of the spacetime manifold undergoes changes, there are too many directions of 
perturbations to include such effects into the initial conditions. Even under this complicated 
situation, we may have several disciplines which make the problem consistent and tractable: 
1. We keep the spherical symmetry; 2. The initial configuration for the scalar fields and the 
corresponding gravitational fields keeps more or less the characteristics of those for n = 1 
static bubble solutions. Furthermore, previous work tells us that the initial size of an 0(4)- 
symmetric bubble at zero temperature is larger than that of an 0(3)-symmetric static bubble 
in the high temperature limit, and that the 0(4)-symmetric bubble expands after a Wick 
rotation. (If a bubble is smaller than the static one, the amount of surface energy required 
to grow is larger than the released bulk energy so that it is energetically favorable to shrink.) 
Thus we assume two initial configurations of the scalar amplitude <fi. One is a scaled n = 1 
bubble solution in which the size of the global monopole is also increased: 

0(O,r) = n=1 (-) (ci>l), (3.1) 

and the other involves the initial shrinking of monopole radius but the expansion of the outer 
bubble wall: 

^»^( [lH.-lAr/^-l)l ) fe>1) ' (3 ' 2) 

Here we solve the equations of motion by numerical calculation. As a coordinate system 
in a Lorentzian spacetime, we adopt the following form, 

ds 2 = -dt 2 + A 2 (t, X )d X 2 + B 2 {t, x)x 2 (d6 2 + sin 2 9d(p 2 ). (3.3) 

With the metric Q3.3| ), we write down the equations of motion ( |2.2j ) and ( |2.3| ) as 

I rrJk P ( N a. 2B ' m 2 ^' -l. 2 ^ -l. dV n (<* a\ 

</, -^-^-(-a + ^ + x^ + 7^ + ¥ = ' (3,4) 
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B' 


(4 


--) 






-4) 


A 2 B 2 




b) 


+ XA 2 


\A 





A 2 B 

x 2 " v ~~ 52 J = (y + 2A 2 ~x 2 b 2 ' (3 - 5) 
1„, /£' 1 



-G* x = ^ + — + -J (3i^ e — K) = AnGH>' (3.6) 

i(G* + G e e + - G\) = K-K 2 + AKK e e - QK 6 * = 8nG(<j) 2 - V), (3.7) 

where the overdot ' and the prime ' stand for the partial derivative with respect to t and x i n 
Eq.( |3.3| ), respectively. Following Ref . |p4fl , we have introduced the extrinsic curvature tensor 



of a t — const hypersurface, whose components are given by 



and we have denoted its trace by K = K\. 

For initial data for the metric, we assume an asymptotically flat de Sitter metric just for 
a technical reason. Although the metric of a Euclidean spacetime is asymptotically closed de 
Sitter spacetime, the effect of the spatial curvature is not so important as long as a bubble is 
smaller than the horizon, and hence our treatment can be verified. Specifically, we suppose 
A(t = 0, r) = B(t = 0, r) = 1 and solve the constraint equations ( |3 .5| ) and ( |3.6| ) to determine 
K{t = 0,r) and K e e {t = 0,r). This treatment is suitable for this system because we obtain 



- f - ~K% - f-f(-T\), (3.9) 

which approaches zero as r increases; we can construct an asymptotically flat spacetime 
without iterative integration. We have also assumed K(t = 0, r) < 0, which means that 
every point in the spacetime is locally expanding. As for the configuration of the scalar field, 
we have supposed Eq.(|3.1|) or ( |3.2| ) and <£>(£ = 0, r) = 0. 

In order to solve the dynamical equations, we use a finite difference method with 1000 
meshes. The evolution of a bubble is depicted by 5 dynamical variables, A, B, K, K e and <j). 
Equations (|3.8|), ( |3.7| ) and ( |3.4| ) provide the next time-step of A, B, K and 0, respectively. 



At each step, we integrate ( |3.6| ) in the radial direction to obtain Kg. In this way we have 
reduced spatial derivatives appearing in the equations, which may become seeds for numerical 
instability. The Hamiltonian constraint equation ( |3.5|) remains unsolved during the evolution 
and is used for checking numerical accuracy. We keep numerical error less than 1 percent 
throughout calculations. 

The results of our numerical computations are summarized as follows. Figure 7 shows that 
the outer wall of an n = 1 bubble starts to expand, whether the bubble wall is thin or thick. 
Figure 7 also tells us that the velocity of the outer wall increases as time elapses. As the 
bubble grows, the wall becomes thinner and the energy accumulated inside the bubble moves 
out to support the expansion of the outer bubble wall (see Fig. 7-(b)). Therefore, after time 
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elapses sufficiently, the motion of the outer wall can be modeled by that of the extremely 
thin wall even for any bubble with initially a thick wall. We trace the time-evolution of 
the position of <fi = 0.5v for the outer wall, as shown in Fig. 8. Its trajectory looks like a 
hyperbola, similar to that of the n = bubble wall. 



(f)/v 




0.5 1 1.5 2 2.5 

R/H- 1 



(f)/v 




0.5 1 1.5 2 2.5 3 

R/H- 1 

Figure 7: n = 1 bubble profiles for fixed A = 1 and v/Mpi = 0.1 associated with the classical 
evolution (ci = 1.2). t/H- 1 = 0, 0.5, 1.0 are shown as dotted, dashed and solid lines, respectively, 
(a) thin-wall bubble of a = 0.1, (b) thick-wall bubble of a = 0.3. 
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Figure 8: Time evolution of the n = 1 bubble radii for A = 1, v/Mpi = 0.1 and c\ = 1.2. (a) 
thin-wall bubble of a = 0.1, (b) thick-wall bubble of a = 0.3. 

An interesting physical quantity at the moment is the terminal velocity in terms of the 
outer expanding coordinates. For an n = bubble it is computed in the thin-wall approxi- 
mation and, because of the gravitational effect, the terminal velocity of the wall is found to 
be smaller than the light velocity |22]. For an n = 1 bubble, because the spacetime between 
the inner wall and the outer wall is described by Eq. fl2.24p and the spacetime outside the 
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outer wall is de Sitter spacetime fl2.25|) , we obtain a junction condition, 



(if) +l 



-nGV(0)R 2 



(if) +1 



8nGv 2 + 



2GIM 



turn | 



R 



-AnGaR, (3.10) 



where R, r, and a are the circumference radius, the proper time, and the surface energy 
density of the shell, respectively e + and e~ are the signs of the angular component of the 
extrinsic curvature of the 2+1-world-hypersurface. In general e + and e~ take +1 or —1, but 
in the present case they are always positive. Neglecting the term of 2GM turn /r in Eq.( |3.10| ), 
we obtain the initial radius of the n — 1 bubble, 



Rn=l(0) 



Rn= (0) 



VI - 8nGv 2 + 



\ 



1 - 8itGv 2 + 



Av 2 



R n =o(0)a 



(3-11) 



where _R n=0 (0) = 3cr/(V(0) + QnGa 2 ) is the initial radius of the n 
of a bubble is sufficiently large, Eq. (|3.11|) reduces to 



Rn=M 



Vl - 8nGv< 



bubble. When the size 



(3.12) 



which shows why the size of R n= i(0) is smaller than i? n= o(0) (Fig. 2-(a)). We numerically 
solve Eq.( p.lO| ) and obtain an approximate formula for the terminal velocity: 



^terminal 



1 - 



4nGV(0) 
3[^o(0)F 



(3.13) 



which agrees with that of an n = bubble. Note that we did not consider the temperature 
effect here in Eq.( |3.13| ), and the inclusion of the temperature effect can decrease the terminal 
velocity further |15[]. Now, let us explain why the long-range term does not change the 
terminal velocity v t0I . minal . When the velocity of the outer wall reaches its terminal velocity, 
an ii = 1 bubble can be approximated as a thin-wall bubble and the energy density T\ 
contributed from the expansion of the long-range tail of the global monopole is 



1 .. 



+ 



2A 2 



+ 



X 2 B< 



+ V 



X 2 B< 



(3.14) 



since the scalar field <fi stays at the vacuum v. This contribution is considerable at the initial 
stage of the evolution; however, it subsides to zero as the n = 1 bubble grows (Bx increases). 
Therefore, the energy difference per unit surface area of the bubble becomes the same as that 
of the n = bubble. We thus conclude that the terminal velocity of an n = 1 bubble is equal 
to that of an n = bubble. 

Another important motion for the n = 1 bubble is of course that of the inner wall. Figure 
8 shows the trajectory of the position of = 0.5v for the inner wall. We find that the inner 
wall just oscillates, and we can expect that this oscillation will be damped gradually. Since we 
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assumed initial configurations with both enlarged (Eq. fl3.1|) ) and shrunken (Eq.( p.2|) ) cores 
of the global monopole, the above result implies that the global monopole inside an n = 1 
bubble is stable against the perturbations of the scalar amplitude, 

f = 0>(x) -> rm, x) + *))• (3.15) 

^From the above discussion, we summarize the evolution of bubbles as follows. When 
the spherical symmetry is assumed for the scalar field, the motion of both n = and n = 1 
bubbles is represented by the expansion of the bubble wall. The global monopole formed in 
the n — 1 bubble remains to be stable and its long-range energy tail in Eq. ( |3.14| ) keeps growing 
before bubble percolation by consuming a part of the false vacuum energy (proportional 
to the increment of bubble radius) obtained from the growth of the true vacuum bubble 
(proportional to the increment of spatial volume). It explains why one need not worry about 
the huge amount of energy necessary to maintain a large-size global monopole if it is formed 
through the first-order phase transition. As a monopole-bubble is larger, the ratio of the 
energy for keeping the global monopole to the energy obtained through the growth of the 
bubble is smaller. Furthermore, this ratio finally becomes negligible for an extremely large 
bubble. We have also shown that the global monopole is a stable object during the evolution 
of a spherically symmetric n — 1 bubble; however, the stability and its physical implication 
due to the distortion of the bubble or the collision of two bubbles remain topics for future 
work. 

One way to understand the stability of the inner bubble wall against the perturbation of 
the scalar amplitude in Eq. ( |3.15| ) is to count the number of negative modes by considering the 
small fluctuation around a given n — 1 bubble solution. Since it is too difficult to calculate 
the negative modes with the inclusion of gravity even for a n = bubble, let us attempt 
to do it in a flat spacetime. Under small static fluctuations around a spherically-symmetric 
bubble solution 8(f)® = </>^ I] CfcV'fc^ 4 )' we obtain a Schrodinger-type equation by varying the 
scalar equation (|2~2"|): 



^2 d 2 V 
-V 2 + — - 



fyktf) = h&M^), (3-16) 



where x % {% = 1,2,3) denote spatial Cartesian coordinates and a subscript n takes or 1, 
corresponding to an n = or n = 1 bubble. For an n = bubble with 0q = (0,0,1), 
Eq.( |3.16|) contains a unique negative mode of which the wave function is that of a nodeless 
s-wave |2|, |l^]. In the case of n = 1 bubbles, Eq. (|3.16|) becomes 



(-v 2 + ^ + ^ )M^) = hM^), (3-17) 

and then the lowest mode is not the nodeless s-wave mode but the I = 1 mode with a single 
node at x = 0. For n = bubbles, the unique negative mode was used for the explanation 
of the motion of its bubble wall. We already showed through the numerical computation of 
n = 1 bubbles that the inner bubble wall is stable but the outer bubble wall starts to evolve 
as time goes. Thus we are likely to interpret this unique negative mode obtained from radial 
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perturbation as the one related with the motion of the outer bubble wall, albeit we need 
further study to reach a definite conclusion. However, for perturbations in all directions, it 
is extremely difficult to count the number of zero modes for an n = 1 bubble in a curved 
spacetime and the n = 1 bubble can also have the possibility of containing more than one 
zero modes as happens for the n = bubble [25 . 



Finally, let us discuss the case where the scale of symmetry breakdown v approaches the 
Planck scale Mp\. When the de Sitter horizon is comparable to or smaller than the radius of a 
bubble, the procedure of a first-order phase transition may not follow the scenario in Ref . || 
but drastic change occurs. For n = bubbles, one may bring up two proposed scenarios 
of phase transitions: One is the one-bubble universe formed inside a thin-wall bubble p3 



and the other is the Hawking- Moss type phase transition p7 |. If we follow the viewpoint 
of Ref. [p3|j for the n = 1 bubble, it is probable that a bubble with a super-horizon-sized 
monopole is nucleated. Once such a configuration is formed, interesting phenomena are 
expected: the evolution of global monopoles [26] or the defect inflation at the monopole site 
0, p. If we consider the model of interest with the 5*0(3) gauge coupling, the nucleated 
false vacuum island is dressed by gauge fields, and consequently shows defect inflation and 
the creation of a Schwarzschild-like wormhole |7|, [TOf. 



IV. Conclusion and Discussions 



In this paper we have studied a first-order phase transition in an 0(3)-symmetric 
model in a curved spacetime and at high temperature. We found a new bubble solution which 
describes another possible decay channel. Different from an ordinary bubble, it contains a 
matter lump at the center of the bubble, which is nothing but a global monopole supported 
by the winding between the internal group 0(3) space and the real space. It manifestly 
shows how the continuous internal symmetry of the theory can play an important role from 
the beginning of bubble nucleation. 

The obtained monopole-bubble (n = 1 bubble) has the following characteristics. First, 
in addition to the outer bubble wall which distinguishes the true vacuum region inside the 
bubble from the false vacuum environment, there is another inner bubble wall which divides 
the core of the global monopole and the true vacuum region inside the bubble with the 
long-range energy tail of the global monopole. Since the formation of the global monopole 
consumes a part of the energy obtained by the difference between false vacuum energy and 
true vacuum energy, the size of the monopole-bubble is slightly larger than that of the 
ordinary bubble for thick-wall bubbles. However, strikingly enough, the opposite is true for 
the monopole-bubbles with sufficiently thin walls. Second, the long-range tail of the global 
monopole is terminated by the cutoff "outer bubble wall" , and the total energy to form the 
global monopole is proportional to the radius of the monopole-bubble. Furthermore, while 
the outer wall expands, the global monopole itself is a stable configuration before bubble 
collisions. It suggests a new mechanism for the production of global monopoles through a 
first-order phase transition despite its infinite energy. The spacetime inside the monopole- 
bubble is flat with solid deficit angle, but the black hole is not produced at the center even at 
the Planck scale. Extending our analysis to general models, we expect that other topological 
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defects or non-topological solitons can also be created through a first-order phase transition as 
solitonic-bubble solutions. We should take into account a new possibility "soliton production 
at the bubble nucleation era" in addition to the soliton production by the horizon or bubble 



collisions 29, 25 



Since the action of a monopole-bubble is larger than that of an ordinary bubble, the 
production rate of monopole-bubbles is exponentially suppressed in comparison with that 
of ordinary bubbles at low temperature; however, it is enhanced considerably at high tem- 
perature. We showed that the production rate of monopole-bubbles can be comparable to 
that of ordinary bubbles for some parameter range of the scalar potential, so the first-order 
phase transition at high temperature described by such models can proceed through this new 
decay channel by the nucleation of monopole-bubbles. If we look at the evolution of the 
monopole-bubble by classical dynamics after the background universe is cooled down to zero 
temperature, we can easily notice the following: 1. The outer bubble wall immediately starts 
to expand and reaches the terminal velocity smaller than light velocity, just as the case of 
an ordinary bubble; 2. The inner bubble wall remains stable and the long-range tail of the 
global monopole grows as the outer wall expands. 

It is well known that topological defects are produced by the Kibble mechanism or bub- 
ble collisions in the early universe [3D], this monopole-bubble nucleation can be a new 
mechanism of producing global monopoles through a first-order phase transition in the early 
universe. If we introduce the gauge coupling, the character of a monopole inside the bubble 
changes from the global one to the 't Hooft-Polyakov monopole [^TJ. It may be interesting 
to apply this mechanism to the production of local monopoles and compare with the results 
obtained in other ways |p5| . If the production rate of monopole-bubbles is too large, then 
one must worry about monopole abundance in the early universe, though it can be diluted by 
inflation. For computing the number of monopoles which survive after the completion of the 
first-order phase transition, further study of bubble collisions is needed, particularly between 
n = 1 and n = bubbles, or between n = 1 and n = 1 bubbles. Finally, we should emphasize 
again that the above procedure to nucleate new bubbles involving a soliton in its center is due 
to continuous symmetry and can be generalized for any continuous global or local symmetry 
which has an appropriate winding between internal group space and spacetime. However, it 
is unclear for the discrete symmetry case in (3+1) dimensions p2| and needs further study. 



Acknowledgments 

The authors are very grateful to G. 't Hooft, V.P. Nair, K. Nakao, M. Peskin, M. Sasaki, 
K. Sato, T. Tanaka, E.J. Weinberg, Piljin Yi, and Zae-young Ghim for valuable discussions. 
Thanks are also due to P. Haines for correcting the manuscript. K.M. and N.S.'s research 
was supported in part by the Grant-in-aid for Scientific Research Fund of the Ministry of 
Education, Science and Culture (No.06302021, No.06640412, and No.07740226), and by the 
Waseda University Grant for Special Research Projects. Y.K.'s research was supported in 
part by JSPS (No.93033), the KOSEF (95-0702-04-01-3, Brain Pool Program) and the Korean 
Ministry of Education (BSRI-95-2413). N.S. thanks Yukawa Institute of Theoretical Physics 
for financial support during his visit. Y.K. thanks Physics Department of Hanyang University 



22 



(BSRI-95-2441) and Center for Theoretical Physics of Seoul National University for their 
hospitality and financial support during his visit. 



References 

[1] D. A. Kirzhnits and A. D. Linde, Phys. Lett. B 42, 471 (1972); L. Dolan and R. Jackiw, 
Phys. Rev. D 9, 3320 (1974); S. Weinberg, ibid 9, 3357 (1974). 

[2] S. Coleman, Phys. Rev. D 15, 2929 (1977); C. Callan and S. Coleman, ibid 16, 1762 
(1977). 

[3] S. Coleman and F. De Luccia, Phys. Rev. D 21, 3305 (1980). 
[4] I. Affleck, Phys. Rev. Lett. 46, 306 (1981). 

[5] A. D. Linde, Phys. Lett. B 70, 306 (1977); ibid 100, 37 (1981); Nucl. Phys. B 216, 421 
(1983). 

[6] A. H. Guth, Phys. Rev. D 23, 347 (1981); K. Sato, Mon. Not. R. Astron. Soc. 195, 
467 (1981); A. D. Linde, Phys. Rev. B 108, 389 (1982); A. Albrecht and P. Steinhardt, 
Phys. Rev. Lett. 48, 1220 (1982). 

[7] K. Sato, M. Sasaki, H. Kodama and K. Maeda, Prog. Theor. Phys. 65, 1143 (1981); K. 
Maeda, K. Sato, M. Sasaki and H. Kodama, Phys. Lett. B 108, 98 (1982). 

A. Vilenkin, Phys. Rev. Lett. 72, 3137 (1994); A. D. Linde, Phys. Lett. B 327, 208 
(1994); A. D. Linde and D. Linde, Phys. Rev. D 50, 2456 (1994). 

[9] N. Sakai, H. Shinkai, T. Tachizawa, and K. Maeda, Phys. Rev. D 53, 655 (1996). 

[10] N. Sakai, Waseda university preprint WU-AP/52/95, |gr-qc/9512045 . 



[11] For a review, see A. Vilenkin, Phys. Rep. 121, 263 (1985) and A. Vilenkin and E. P. 
S. Shellard, Cosmic Strings and Other Topological Defects, (Cambridge, Cambridge 
University Press, 1994). 

[12] S. Coleman, Nucl. Phys. B 298, 178 (1988). 



[13] Y. Kim, Nagoya university preprint DPNU-94-39, |hep-th/9410"076 . 



[14] M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63 341 (1989); D. Harari and C. Lousto, 
Phys. Rev. D 42, 2626 (1990). 

[15] P. J. Steinhardt, Phys. Rev. D 25, 2074 (1982). 

[16] G. H. Derrick, J. Math. Phys. 5, 1252 (1964); R. Hobart, Proc. Phys. Soc. 82, 201 
(1963). 



23 



R. Gregory, Phys. Lett. B 215, 663 (1988); A. G. Cohen and D. B. Kaplan, ibid 215, 
67 (1988); D. Harrari and P. Sikivie, Phys. Rev. D 37, 3438 (1988). 

D. Harrari and A. Polychronakos, Phys. Lett. B 240, 55 (1990); G. W. Gibbons, M. E. 
Ortiz and F. Ruiz Ruiz, ibid B 240, 50 (1990). 

For the rigorous proof of existence of bubble solutions in flat spacetime, see Ref. || for 



an n = bubble, and Ref. [13] with the help of C. Kim, S. Kim and Y. Kim, Phys. Rev. 



D 47 (1993) 5434 for an n = 1 bubble. 

J. R. Gott, Astrophys. J. 288, 422 (1985); B. Linet, Class. Quantum Grav. 7, L75 
(1990); M. E. Ortiz, Phys. Rev. D 43, 2521 (1991); C. Kim and Y. Kim, ibid 50, 1040 

(1994) . 

C. J. Hogan, Phys. Lett. B 133, 172 (1983); M. S. Turner, E. J. Weinberg and L. M. 
Widrow, Phys. Rev. D 46, 2384 (1992). 

V. A. Berezin, V. A. Kuzmin and I. I. Tkachev, Phys. Lett. B 120, 91 (1983), Phys. 
Rev. D 36, 2919 (1987); S. K. Blau, E. I. Guendelman and A. H. Guth, ibid 35, 1747 
(1987); N. Sakai and K. Maeda, ibid 50, 5425 (1994). 

A. H. Guth and E. J. Weinberg, Nucl. Phys. B 212, 321 (1983). 

R. Arnowitt, S. Deser and C. W. Misner, in Gravitation: An Introduction to Current 
Research edited by L. Witten, (Wiley, New York, 1962). 

R. Basu, A. H. Guth and A. Vilenkin, Phys. Rev. D 44, 340 (1991); J. Garriga and A. 
Vilenkin, ibid 47, 3265 (1993). 

E. I. Guendelman and A. Rabinowitz, Phys. Rev. D 44, 3152 (1991). 

S. W. Hawking and I. G. Moss, Phys. Lett. B 110, 35 (1982); L. G. Jensen and P. J. 
Steinhardt, Nucl. Phys. B 237, 176 (1984), ibid B 317, 693 (1989); D. A. Samuel and 
W. A. Hiscock, Phys. Rev. D 44, 3052 (1991). 

K. Lee, V. P. Nair and E. J. Weinberg, Phys. Rev. Lett. 68, 1100 (1992), Phys. Rev. 
D 45, 2751 (1992); M. E. Ortiz, ibid 45, R2586 (1992); P. Breitenlohner, P. Forgacs, 
and D. Maison, Nucl. Phys. B 383, 357 (1992); T. Tachizawa, K. Maeda and T. Torii, 
Phys. Rev. D 51, 4054 (1995). 

T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 

J. Borrill, T. W. B. Kibble, T. Vachaspati and A. Vilenkin, Phys. Rev. D 52, 1934 

(1995) . 

Y. Kim, K. Maeda and N. Sakai, in preparation. 

I. V. Barashenkov and E. Yu. Panova, Physica D 69, 114 (1993). 



24 



